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Введение. Статья посвящена исследованию работы гидро- 
динамического кавитатора, используемого при эрозийном 
воздействии на поверхность твердого тела, а также опти- 
мизации структуры устройства для увеличения разруши- 
тельной способности кавитационной струи. В работе рас- 
смотрено влияние единичного дефекта внутренней по- 
верхности комбинированного сопла на объемную долю 
содержания пара и геометрию области кавитации. Целью 
работы является выявление методом численного модели- 
рования закономерностей влияния дефекта внутренней 
поверхности кавитатора различной величины на гидроди- 
намические и кавитационные характеристики сопла. 
Материалы и методы. Использованы возможности про- 
граммного пакета конечно-элементного анализа АМЗУ5 
У/огКБепсв и интегрированного в него модуля оптимиза- 
ции процесса разработки и технологической подготовки в 
области вычислительной динамики жидкостей и газов 
АМЗУЗ СЕХ. В основу моделирования положены экспе- 
риментальные данные, полученные при истечении воды в 
кавитационном режиме из исследуемого сопла на специ- 
ально разработанном лабораторном стенде. 

Результаты исследования. Получены и представлены 
графические зависимости объемной доли содержания 
пара, полного давления и длины кавитационной области 
от расстояния вдоль оси струи при различной величине 
дефекта. Выявлены две фазы течения кавитационной 
струи в неидеальном кавитаторе, и показано влияние этого 
перехода на распределение скоростей в сечении устрой- 
ства. 

Обсуждение и заключения. Наличие внутреннего дефекта 
на поверхности расходящегося конусного участка комби- 
нированного сопла размером менее четверти диаметра 
центрального цилиндрического участка может не вызы- 
вать визуальных изменений в геометрии области кавита- 
ции, однако значительно снижает эрозийную способность 
кавитационной струи. Дальнейшее увеличения дефекта 
приводит к полному подавлению кавитации потока, но 
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Введение. Кавитация представляет собой многофазное и сложное физическое явление, определяемое 
как образование пузырьков пара в жидкости, когда давление локально падает ниже давления насыщенного па- 
ра. Рост пузырьков пара в метастабильной жидкости был впервые исследован Рэлеем [1], который представил 
приближенное аналитическое решение для роста сферически-симметричного пузырька только с учетом инер- 
ционной силы. Впоследствии исследования по процессу роста пузырьков проводились и проводятся во многих 
работах [2-4]. 

Кавитация классифицируется безразмерным параметром — числом кавитации: 

2(р-р,) 
ро’ 
гдер — противодавление в затопленной полости, р, — давление насыщенных паров жидкости при температуре 


окружающей среды, р — плотность жидкости, > — скорость потока. С уменьшением числа кавитации, путем 
увеличения скорости жидкости, наблюдается суперкавитация, которая заключается в образовании устойчивой 


У = 


газо-паровой области. 

В ранние годы кавитация была нежелательным явлением, поскольку наносила большой ущерб гидрав- 
лическим машинам, турбинам и насосам, что приводило к разным уровням потери их эффективности в зависи- 
мости от интенсивности кавитации. Однако с увеличением числа исследований по кавитационному поведению 
было доказано, что кавитационные струи могут быть практически применимы во множестве полезных прило- 
жений. Кроме уже широко используемого явления ускорения химических реакций при приготовлении эмуль- 
сий и суспензий [5, 6], появляются новые современные технические решения использования суперкавитации. 
Среди них широко исследуемая экспериментально и теоретически [7, 8], в том числе и с использованием про- 
граммных пакетов моделирования гидродинамических процессов [9-12], проблема снижения сопротивления 
движению подводных высокоскоростных аппаратов. Отдельно следует выделить область разработки и проек- 
тирования дизельных тепловых двигателей, в которых кавитация применяется не только для приготовления 
водно-дизельной эмульсии [13], но и серьезно влияет на коэффициент расхода топлива систем впрыска под вы- 
соким давлением [14]. 

Применение кавитационных струй главным образом связано с локальным, резким увеличением давле- 
ния и температуры, возникающим в результате коллапса кавитационных пузырьков, а также с довольно силь- 
ными колебаниями давления, вызванными соударением кавитационной струи с твердой поверхностью. В более 
конкретных выражениях может быть получено давление в приближении 10 ГПа и температура выше 5000 К во 
время пикосекундного коллапса кавитационных пузырьков [15], что вызывает быстрое разрушение хрупких 
материалов. Более того, процесс схлопывания кавитационных пузырьков способен возбуждать ударные волны 
и высокоскоростные микроструйки, которые считаются основной механической причиной сильной эрозионной 
способности кавитационных струй [16-18]. По этой причине кавитационную эрозию применяют при очистке 
подводных конструкций и механизмов от загрязнений, биологических обрастаний и прочих нежелательных 
образований [18-22]. При эксплуатации морских судов, базирующихся в теплых морях, возникает большая 
проблема обрастания их подводной части различными моллюсками и балянусом. Наиболее эффективным и со- 
временным способом очистки от обрастаний поверхностей под водой является гидродинамический, использую- 
щий как динамическое воздействие, так и кавитационный эффект. В некоторых случаях, например, при обработ- 
ке пластмассовых и стеклопластиковых корпусов, кавитационная очистка — единственное возможное решение 
во избежание повреждений поверхности. 
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Методы усиления интенсивности кавитации струй привлекли к себе значительное внимание в течение 
последних нескольких лет. При решении этого вопроса многие исследователи фокусируются на механизмах 
генерации, роста и коллапса кавитационных пузырьков, геометрии кавитационной области за счет оптимизации 
структуры сопел для создания более разрушительных кавитационных струй [23]. 

Для изготовления сопловых насадок и кавитаторов, которые в дальнейшем используются для создания 
струйных потоков жидкости, истекающей из них, требуется специальное оборудование. Это должны быть пре- 
цизионный станок, способный изготовить изделие с наименьшим биением шпинделя, специальный инструмент 
для данного вида сопловой насадки и строгое соблюдение технологического процесса, разработанного механи- 
ком-технологом. Чаще всего при изготовлении сопловых насадок не соблюдается методика их изготовления, и 
они не отвечают требованиям технического задания. 

Долгое время исследование влияния дефектов и шероховатости внутренней стенки кавитаторов на по- 
ток и кавитационную характеристику сопла рассматривались как бесполезный параметр в общих исследовани- 
ях. Поэтому, несмотря на большое количество работ выполненных по моделированию геометрий или структур 
сопел, исследование дефектов обработки внутренней поверхности сопла и их влияние на кавитационною эро- 
зию еще не было выполнено. Вместе с тем, шероховатость поверхности может иметь значительные последствия 
для генерации кавитационных пузырьков, длины кавитационной области, распределения объемной доли пара. 

Кроме того, исследования кавитации в сопле инжектора дизельного двигателя показали, что при увели- 
чении шероховатости внутренней поверхности, массовый расход и коэффициент расхода сначала возрастают, а 
затем уменьшаются [14]. Основываясь на анализе вариаций объемной доли пара и средней скорости потока на 
выходе из отверстия, было сделано заключение, что определенный уровень черновой обработки внутренней 
стенки отверстия может изменить характеристики потока за счет значительного подавления кавитации внутри 
сопла. 

Таким образом, цель данной работы — исследовать методом численного моделирования влияние нали- 
чия дефектов внутренней поверхности кавитатора на кавитационные характеристики сопла, используемого для 
эрозийного воздействия при обработке твердых материалов. 

Материалы и методы 

Экспериментальные данные. С помощью визуальных наблюдений истечения струйных потоков из 
сопловых насадок можно частично выявлять некачественно изготовленные изделия. Однако окончательную 
оценку качественного изготовления можно дать только после комплекса испытаний изделия на специально раз- 
работанном стенде [24, 25], оборудованном чувствительными датчиками и контрольной аппаратурой. 

Сопловые насадки, выполняющие роль возбудителей кавитации (кавитаторы), наиболее чувствительны 
к качеству их изготовления. При различных значениях гидростатического давления в окружающей затоплен- 
ную струю жидкости кавитационный поток, истекающий из кавитатора изготовленного некачественно, может 
изменять свою структуру, направление движения, что существенно снижает эрозионную способность кавита- 
ЦИИ. 

На рис. 1 представлены фотографии затопленных струй, истекающих из кавитатора, сделанного из 
прозрачного материала. Давление, подводимое к кавитатору — р»=1О0МПа, а противодавление в затопленной 
около кавитатора области изменяется в интервале р,=0,1-0,4 МПа. На рис.1 хорошо просматривается тот факт, 
что с увеличением противодавления происходит смещение зоны скопления газопаровых пузырьков (каверны) 
струйного потока вправо. Этот факт позволяет судить о некачественном изготовлении внутренней поверхности 
кавитатора. Разработанный гидрокавитационный комплекс [24, 25] позволяет обнаружить следствие наличия 
дефекта обработки сопловой насадки. Причины искажения формы кавитационной струи можно выяснить при 
рассмотрении заведомо известной погрешности изготовления. Однако точно задать геометрию, размер и 
расположение нарушения плоскостности внутренней поверхности кавитатора — не простая практическая 
задача. Решить подобную проблему позволяет исследование гидродинамической кавитации при помощи 


компьютерного моделирования струйного истечения. 
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©) 
Рис. 1. Вид на истечение затопленных струй при различных параметрах 
противодавления в окружающей струю жидкости: 


а — противодавление рх = 0,1 МПа, 6 — р = 0,2МПа, с — рк= 0,4 МПа 


Вр. 1. Иез оГиптегуе4 уе @эсйагее ай уапоиз Баску’агА ргеззиге рагатеетх 
т атЫет! Паша: а —фБасКи’агА ргеззиге ру = 0.1МРа, Ь — ру = 0.2МРа, с — рь= 0.4МРа 


Моделирование гидродинамических процессов средствами пакета конечных элементов АМЗУЗ$ СЕХ — 
инструмент, который уже широко используется для оптимизации и проектирования различных частей гидроди- 
намических устройств и их многофункционального анализа. АМЗУ$ СЕХ является высокопроизводительной 
программой вычислительной гидродинамики, обеспечивающей надёжное и быстрое решение широкого круга 
задач, связанных с течениями жидкостей и газов, в том числе и с учетом кавитации [26-29]. 

Методика численного моделирования. При моделировании использована функция кавитации в АМЗУ$ 
СЕХ, которая основана на уравнении Рэлея-Плессета: 

2 -2° _Рн-Р 
АР З\ 4% рК р 
где К — радиус кавитационного пузырька, р, — давление внутри пузырька (в модели — давление насыщенных 
паров), р — локальное давление в жидкости (абсолютное давление СЕ)-решателя), р — плотность жидкости, с 
— коэффициент поверхностного натяжения между жидкостью и паром. 

Поскольку решить уравнение (1) в общем виде невозможно, в модели Рэлея-Плессета АМЗУЗ СЕХ его 
упрощают, пренебрегая вторым порядком и не учитывая поверхностное натяжение жидкости: 





В (1) 





4 _ |[2Рн-Р 
Ги З р ‘ 
Скорость изменения объема пузырьков пара рассчитывается как 
ат _а 4 пр — 4тВ? ай. 
4 3 3 р 


а скорость изменения массы 





Уравнения, выражающие соотношения масс выделившегося и сконденсированного пара, принятые в 
расчете, имеют вид [10]: 











Зар, |2 рн-Р 
т =Е* & Н ( < ) 
в р (Р<Ры 
-_ 30% (1-@)р, [2Ри-Р 
ай - - > 
т г т (Р> ры), 


где и — коэффициент связи объемной доли пара с массовой; К — начальный радиус пузырька, 
р, — плотность паровой фазы, Е", Е` — константы моделирования. 

Модель Рэлея-Плессета учитывает два фактора: скорость роста пузырька и статистический характер 
распределения пузырьков в кавитационном потоке. Объемная доля пара в ячейке учитывается статистически на 
основе определения числа зародышей п и их радиуса Ау и определяется по формуле: 


И - 
а=-лАп. 
3 КЮ 
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Модель кавитации Рэлея-Плессета, реализуемая в АМЗУЗ СЕХ, требует задания таких параметров, как 
средний диаметр пузырька и давление насыщенного пара. По данным разных авторов [30-32], радиус пузырька 
в конце сжатия Ани» как правило, находится в пределах 10—10 м при начальном радиусе А\=1х10°м. Таким 
образом, объем пузырька изменяется не менее, чем в 1000 раз, что и вызывает экстремальную концентрацию 
первоначально запасенной пузырьком энергии. Основной диаметр пузырьков пара (Маш Фатеег) в уравнении 
Рэлея-Плессета выбран 2х 10° м, давление насыщенных паров (ЗагаНоп Ртеззиге) с учетом температуры моде- 
лируемой жидкости 25°С р„ = 3170 Па. Остальные параметры моделирования, выбранные по умолчанию, счи- 
таются стандартными для большинства кавитационных течений. 

В качестве исходного образца выбрано комбинированное сопло, состоящее из трех участков (рис. 2а). 
Первый — конусный сходящийся участок с углом 01=14° заканчивается цилиндрической частью внутренним 
диаметром 4 = 2мм и длиной 5 = 1,754. Третий участок конусной формы имеет угол раскрытия 9›=45°. Истече- 
ние жидкости происходит в затопленное пространство (рис. 2 5), противодавление в котором задается в пре- 
процессоре СЕХ и для всех расчетных моделей составляет 0,5 МПа. Сеточная модель тетраэдрического типа 
сгенерирована общим числом элементов 1,28 млн. Исследуемая жидкость — вода, плотностью р=1000 кг/м* 
при температуре 7=25°С. Массовый расход через гидродинамическую систему 0,45 кг/с выбирался с учетом 
того, что полное давление по оси струи при реальном эрозийном использовании кавитационного потока состав- 
ляет р, = 10-15 МПа. 





0,00 40,00 (ппгп) 
20,00 





а) Ь) с) 
Рис. 2. Геометрия комбинированного сопла 
а — схема проточной части кавитатора; Ь — внешний вид модели; 
с — внутренний дефект расходящегося конусного участка 
Ею. 2. СотЬте4 по22 эеотегу 
а — Чаэтат офсауйиог Паш епа; Ь — то4е! арреагапсе; 
с — пиегпа! 4еесй оГсотЬте4 пог сотса[ райетп 


Наличие повреждений моделировалось на внутренней поверхности расходящегося конусного участка 
(рис. 2 с). Дефект в сечении представлял равнобедренный треугольник с углом 0;=45° между равными сторона- 
ми 6. Вершина треугольника с углом 6 может находиться снаружи или внутри выходного конуса, тем самым 
уменьшая его объем. 

Результаты исследования. На рис. 3 а представлен график зависимости объемной доли содержания 
пара от расстояния вдоль оси струи х для различных значений 6 стороны треугольного дефекта. Начало коор- 
динат соответствует началу цилиндрического участка. Для всех д зависимость и = Кх) имеет максимум, кото- 
рый с увеличением 6 приближается к цилиндрическому участку. В случае с идеальным соплом (кривая 1 на 
рис. 3 а) максимум содержания паровой фазы в объеме приходится на значение и = 0,6, которое незначительно 
меняется при наличии в выходном конусе наружного дефекта (кривая 2 на рис. За). С увеличением объема рас- 
ходящегося конуса за счет треугольного выступа распределение у = Хх) заметно отличается от идеального слу- 
чая только во второй части кавитационного «факела». Это отклонение в основном связано с изменением формы 
и не влияет на максимум концентрации пара и длину кавитационной области. Поэтому на рис. 3 а показана 
только одна зависимость и’ = Кх) сопла с наружным дефектом 6=0,5мм. Для остальных значений 6 кривые на 
графике плотно располагаются и мало отличимы. 
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Рис. 3. Зависимость объемной доли содержания пара у» (а) и полного давления р, (Ъ) от расстояния вдоль оси струи х при 





различной величине дефекта 6, измеряемого в мм: 1 — 0; 2 — 0,5 (дефект наружу); 3 — 0,3; 4 — 0,4; 5 — 0,5; 
6 — 0,7; 7 — 0,8; 8 — 0,85;9 — 0,9; 10 —1 





Е1о. 3. ереп4епсе оР уоте йасйоп оЁ уарог сощепе и’ (а) ап 1ю1а ргеззиге р, (Ъ) оп Чз$апсе аопо ]е! Х-ах1$ Юг уап!- 
оиз деРесЕ уаез д ш тит; 1-0; 2—0.5 (4еесЁ об); 3—0.3; 4—0.4; 5-0.5; 6-0.7; 7—0.8; 8—0.85; 9—0.9; 10-1. 





Образование треугольного дефекта с вершиной внутрь объема расходящегося конуса снижает развитие 
кавитации, а при величине б = 0,54 полностью подавляет ее. Максимум содержания пара в потоке жидкости 
падает и при 6 = 1мм составляет менее 0,001%. С увеличением д длина кавитационной зоны / меняется нерав- 
номерно (рис. 4) и по характеру спада не повторяет у. При значениях д < 0,24 [ незначительно возрастает, а в 
интервале 0,2 < 6 < 0,44% остается равной длине кавитационного «факела» идеального сопла. Дальнейшее уве- 
личение 6 приводит к резкому сокращению области кавитации и [=0 при 6=0,54%. 
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Рис. 4. Зависимость длины кавитационной области / от величины дефекта 6 
Ев.4 Оерепаепсе оГсауйайоп гезоп [епай [ оп аеесЕ уаие д 


Подобная закономерность наблюдается и для зависимости полного давления потока жидкости р, от 
расстояния вдоль оси струи (рис. 3 Б). Поскольку постоянной величиной в моделировании принят массовый 
расход, полное давление является функцией переменной, зависящей от кавитации струйного истечения. Нали- 
чие дефектов в выходном конусе сопла размером д < 0,44 не влияет на величину р. на цилиндрическом участке 
(кривые 1-7 на рис. 3 Ь) и качественно повторяет зависимость р. = Дх) идеально выполненного кавитатора, имея 
резкий спад после начала зарождения кавитации. Он достигает своего минимума в точке наибольшего значения 
и’ (рис. 3 а). Далее происходит плавное выравнивание р, независимо от величины дефекта. Если же д > 0,44), 
входное давление заметно возрастает (кривые 8—10 на рис. 3 Ь) и на протяжении всего сопла остается постоян- 
НЫМ. 

Таким образом, можно выделить две фазы течения кавитационной струи в неидеальном кавитаторе. На 
первом этапе, когда величина внутреннего дефекта менее 0,44, активное зарождение кавитационных пузырь- 
ков начинается непосредственно в конце цилиндрической части сопла, а длина кавитационной области слабо 
зависит от значения 5. Появление выступа, по сравнению с идеальным конусным участком, снижает объемную 
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долю содержания пара и спад полного давления в области зарождения кавитации. Второй этап течения при 
наличии дефекта д > 0,44 характеризуется быстрым уменьшением длины области кавитации и концентрации 
паровой фазы практически до нуля. Уже при б = 0,54 кавитация полностью подавляется внутри сопла, при этом 
полное давление возрастает внутри кавитатора и на выходе из него. 

Кавитация в потоке жидкости наблюдается в местах быстрого изменения скорости течения жидкости, 
которое влечет за собой уменьшение давления. При постоянных значениях массового расхода величины проти- 
водавления в камере истечения струи для несжимаемой жидкости распределение скоростей потока будет опре- 
делять характер зарождения, роста и схлопывания кавитационных пузырьков. На рис. 5 представлены поля ско- 
ростей жидкости для фазы течения с развитой кавитационной областью в идеальном кавитаторе (рис. 5 а) и при 
наличии выступа размером 6 = 0,254 (рис. 5 Б), а также для переходного процесса 6 = 0,44 (рис. 5с) и подавле- 
ния кавитации при 6 = 0,54 (рис. 5 4). Расположение области максимальных скоростей изменяется от величины 
дефекта кавитатора, при этом остается постоянным на уровне 200-180 м/с. Для случая развитой кавитации 
жидкость достигает максимальной скорости в конце цилиндрического участка и уменьшается в расширяющем- 
ся участке гидродинамического устройства. Появление препятствия потоку приводит к появлению второй об- 
ласти повышенных скоростей, которая с увеличением размеров дефекта преобладает над снижающейся скоро- 
стью на входе в конусный участок. 

Согласно уравнению (1) для несжимаемой жидкости процесс кавитации определяется локальным дав- 
лением жидкости, которое в СЕ)-решателе является абсолютным и равно 


Рав > Руци Рег | (2) 


где ри, — статическое давление жидкости, р‚.г—— опорное давление, относительно которого измеряются любые 
давления, установленные в моделировании. В свою очередь, полное давление равно 


р ЯРуки * Рауп’ 2 


2 
где р в ть. динамическое давление жидкости, которое при постоянных значениях максимальной скоро- 
сти о не вносит вклад в падение полного давления для дефекта 6 < 0,44, (кривые 1—6 на рис. 3,6). Распростране- 
ние струи в затопленном пространстве при зарождении кавитации происходит в среде с переменной плотно- 
стью, зависящей от концентрации пузырьков газа или пара. При понижении статического давления на выходе 
из цилиндрического участка кавитатора до величины давления насыщенного пара р, резко увеличивается объем 
газо-жидкостной смеси за счет образования пузырьков, что в изотермических условиях течения приводит к до- 
полнительному снижению давления и усилению кавитации. Появление препятствия потоку уменьшает скорость 
течения в цилиндрическом участке и исключает возможность появления кавитации, нивелируя и падение дав- 
ления при образовании пузырьков пара, что согласно (2) и (3) отображается на зависимостях и = Дх) и р. = Хх) 


(рис. 3). 
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а Рис. 5. Распределение скорости потока жидкости в среднем сечении кавитатора при различных значениях стороны 


треугольного дефекта: а —6=0, 6 — д = 0,5 мм, с —д= 0,8 мм, 4—9=1 мм 


Ер. 5. Глятфиноп оГ ша Лоу? гие т сауйшог табеспоп а! уапои5 ущиез оГтапешаг 4еесй 54е: а- 9 =0,Б-д= 0.5 тт, 
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Обсуждение и заключение. Влияние единичного дефекта на внутренней поверхности расходящегося 
конусного участка на гидродинамические параметры истечения струи воды в кавитационном режиме можно 
разделить на две стадии. При небольших размерах (менее 0,4 части от диаметра цилиндрического участка кави- 
татора) сохраняются геометрические параметры кавитационной области, однако концентрация паровой фазы и 
полное давление уменьшаются. Наличие таких дефектов трудно выявить визуально без дополнительных коли- 
чественных измерений. Они не нарушают симметрию струи, но снижают ее эрозийную способность. Во втором 
случае, если внутренний выступ относительно большой (более 0,4 части от диаметра цилиндрического участка 
кавитатора), происходит резкий спад объемной доли паровой фазы в затопленной струе до полного прекраще- 
ния кавитации. При этом на выходе из гидродинамического устройства полное давление струи возрастает ввиду 
отсутствия снижения статического давления за счет образования пузырьков пара. 

Таким образом, полученные в работе результаты способствуют усовершенствованию конструкций гид- 
родинамических кавитаторов и улучшению их эрозийного воздействия при использовании кавитации для 
очистки подводных конструкций и механизмов. 
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